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i ' Abstract 

A search for squarks of -R-parity violating supersymmetry is performed in ep colli- 
sions at HERA using HI 1994 e + data. Direct single production of squarks of each 
generation by e + -quark fusion via a Yukawa coupling A' is considered. All possible 
^ ■ i?-parity violating decays and gauge decays of the squarks are taken into account. 

No significant deviation from the Standard Model predictions is found in the various 
multi-lepton and multi-jet final states studied and exclusion limits are derived. At 
95% confidence level, the existence of first generation squarks is excluded for masses 
up to 240 GeV for coupling values A' > ^/47ra em . The limits obtained are shown to 
be only weakly dependent on the free parameters of the Minimal Supersymmetric 
Standard Model. Stop squarks are excluded for masses up to 138 GeV for coupling 
A' x cos#t to e + d pairs > 0.1 x v / 47ra eTO , where 6t is the mass mixing angle. 

Light stop squarks are furthermore searched for through pair production in 7- 
gluon fusion processes. No signal is observed and exclusion limits are derived. 
Masses in the range 9 to 24.4 GeV are excluded at 95% confidence level for A' x 
cos 9 t > 10" 4 . 
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1 Introduction 



The search for squarks, the scalar supersymmetric (SUSY) partners of the quarks, is 
especially promising at the ep collider HERA if they possess a lepton number violating 
Yukawa coupling A' to lepton-quark pairs. Such squarks, present in the i?-parity violating 
(ftp) SUSY extension of the Standard Model (SM), can be singly produced via the coupling 
A' as s-channel resonances. Masses up to the kinematic limit of ^/s ~ 300 GeV are 
accessible by the fusion of the 27.5 GeV initial state positron with a quark of the 820 GeV 
incoming proton. In the low mass range, pair production via 7-gluon fusion provides a 
complementary search largely insensitive to the Yukawa coupling. 

In this paper, squarks are searched through single production via a fi p coupling, con- 
sidering both decays and all possible decays via gauge couplings involving mixed states 
of gauginos and higgsinos. A search for pair production of light stops at low masses via 
7-gluon fusion is also carried out. The analysis uses the 1994 e + p data corresponding to 
an integrated luminosity of Cdata = 2.83 pb . Earlier squark searches at HERA were 
presented in |]J. 



2 Phenomenology 

The general SUSY superpotential allows for gauge invariant terms with Yukawa couplings 
between the scalar squarks (q) or sleptons (I) and the known SM fermions. Such couplings 
exist if one assumes the possibility of violating (multiplicatively) the conservation of im- 
parity which is imposed in the Minimal Supersymmetric Standard Model (MSSM); R p = 
(— 1) 3B+L+25 where S denotes the spin, B the baryon number and L the lepton number 
of the particles. Of particular interest for HERA are the $ p terms X^^LiQjD^ of the 
superpotential which allow for lepton number violating processes. By convention the 
ijk indices correspond to the generations of the superfields Lj, Qj and Dk containing 
respectively the left-handed lepton doublet, quark doublet and the right handed quark 
singlet. Expanded in terms of matter fields, the interaction Lagrangian reads Pfl : 



^-LiQjDk — Kjk 



-eX4 - ei«14 - (4)X4* 

* 1 M Ak 1 f-i \cjj jk* 



+P L did% + v L d\d k R + (^M4*J + h-c. 

where the superscripts c denote the charge conjugate spinors and the * the complex con- 
jugate of scalar fields. For the scalars the 'R' and 'L' indices distinguish independent 
fields describing superpartners of right- and left-handed fermions. Hence, with an e + 
in the initial state, the couplings A' ljfc allow for resonant production of squarks through 
positron-quark fusion. The list of possible single production processes is given in table [1|. 
In this paper, the squark search is carried with the simplifying assumptions that: 

• only one of the A' ljA , dominates; 

• squarks and q^) of the first and second generation are quasi-degenerate in mass 
(the case of the stop squark is considered separately) ; 

• the lightest supersymmetric particle is the lightest neutralino x?; 
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Table 1: Squark production pro- 
cesses at HERA (e + beam) via a im- 
parity violating X' x - k coupling. 



ijk 



111 
112 
113 
121 
122 
123 
131 
132 
133 



production process 



+ u 
+ u 
+ u 
+ c- 
+ c- 
+ c 

+ t- 



d R 
• §r 
■Ir 
d R 
sr 
b_R 
d~R 
s R 
b R 



e + + d ■ 
e + + s - 
e + + b- 
e+ + d 
e+ + s 
e+ + b - 
e+ + d 
e+ + s 



u l 
u l 
■ c L 

CL 

cl 

■t L 
■t L 

■h 



• gluinos are heavier than the squarks such that decays q — > q + g are kinematically 
forbidden. 

The squarks decay either via their Yukawa coupling into fermions, or via their gauge 
couplings into a quark and either a neutralino Xi {i — 1)4) or a chargino Xj (j — 1)2). 
The mass eigenstates Xi an d Xj are mixed states of gauginos and higgsinos and are in 
general unstable. In contrast to the MSSM, this also holds in fi p SUSY for the lightest 
super symmetric particle (LSP) which decays via A' ljfe into a quark, an antiquark and a 
lepton g. 

Typical diagrams for the production of first generation squarks are shown in Fig. [1|. 
By gauge symmetry only the d R and ul are produced via the A' couplings. These have in 
general widely different allowed or dominant decay modes. 

In cases where both production and decay occur through a X[ fc coupling (e.g. Fig. [IJa 

and c for A' m ^ 0), the squarks behave as scalar leptoquarks M, |J. For A' m ^ 0, the d R 
resemble the S° leptoquark and decays in either e + + u or u e + d while the ul resemble the 
5*1/2 and only decays into e + d. Hence, the final state signatures consist of a lepton and 
a jet and are, event-by-event, indistinguishable from the SM neutral (NC) and charged 
current (CC) deep inelastic scattering (DIS). The strategy is then to look for resonances 
in DIS-like events at high mass, exploiting the characteristic angular distribution of the 
decay products expected for a scalar particle. 

In cases where the squark decay occurs through gauge couplings (e.g. Fig. |l|b and 
d), one has to consider for the u L the processes Ul — > ux® or dxf while for the d R only 
d R — > o?x° is allowed. This is because the SU(2) l symmetry which implies in the SM that 
the right handed fermions do not couple to the W boson also forbids a coupling of d R to 
the W. Hence, the d R can only weakly couple (in proportion to the d quark mass) to the 
Xj through its higgsino component. 

The possible decay modes of the chargino, when it is the lightest chargino Xii are 
the gauge decays xt ~> Xi^ +hl an d Xi XiQQ'j an d the $ p decays xl ~^ pv >d and 
Xi — > e + dd. The fate of the Xi depends on its gaugino-higgsino composition. The 
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u, d 



f ■ 



— \ 

f \ 



< 



(c) 



(d) 



Figure 1: Lowest order s-channel diagrams for first generation squark production at HERA 
followed by (a),(c) Ijt p decays and (b),(d) gauge decays. In (b) and (d), the emerging neutralino 
or chargino might subsequently undergo IjL v decays of which examples are shown in the doted 
boxes for (b) the \\ an d (d) the Xi- 

question of how this Xi nature depends on free fundamental parameters of the MSSM, as 
well as the corresponding q branching fractions for various possible decay channels will 
be discussed briefly in relation to our analysis in section |5] and was studied in more detail 
in || ||. In general, the Xi wm undergo the decay Xi ~~ > e ± qq' or Xi ~- > vqq. The former 
will be dominant if the Xi is photino-like (i.e. dominated by photino components) in 
which case both the "right" and the "wrong" sign lepton (compared to incident beam) 
are equally probable leading to largely background free striking signatures for lepton 
number violation. The latter will dominate if the Xi is zino-like. A higgsino-like Xi wu l 
most probably be long lived and escape detection since its coupling to fermion-sfermion 
pairs (e.g. Fig. |l|d) is proportional to the fermion mass [0. Hence processes involving a 
if-like Xi wn l De affected by an imbalance in transverse momenta. 

Taking into account the dependence on the nature of the x?, the possible decay chains 
of the ul and cIr squarks can be classified into eight distinguishable event topologies listed 
in tables and |3] and labelled SI to S8. 
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Channel 



xi 

nature 



Decay processes 



Signature 



SI 



1,Z,H 



Q 



High P T e+ + 1 jet 



S2 



1,Z,H 
H 



d R 
Q 





d 






Q 


xi 






Q 


X<t 








X' 

c — > 


e 


+ q'q" 


d 


xt 








X' 


e 


+ dd 


d 


xt 










W+xi 








X' 

c — ^ 








^ q q 



Missing Pj- + 1 jet 



S3 



7^ 



9 



High P T e+ 
+ multiple jets 



e + q'q" 



S4 



1,Z 
1,Z 



Q Xi 



d Xi 



e q'q" 



W + X °i 

A' 



e q'q" 



High Pt e~ 

(i.e. wrong sign lepton) 
+ multiple jets 



q q' 



Table 2: Squark decay channels in Ijt p SUSY classified per distinguishable event topologies 
(first part). The dominant component of the Xi f° r which a given decay chain is relevant is 
given in the second column. The list of processes contributing to a given event topology is 
here representative but not exhaustive, e.g. the gauge decays of the xf involving a virtual W + 
(Fig. ||b) may also proceed via a virtual sfermion. 



For a squark decaying into a quark and the lightest neutralino, the partial width can 
be written as 



k( A2 + B2 ) M *( ] 



M , M\Y 2e 2 el 



q 




where A and B in the left expression are chiral couplings depending on the mixing param- 
eters. Detailed expressions for such couplings can be found in 0. Under the simplifying 
assumption that the neutralino is a pure photino 7, this gauge decay width reduces to 
the expression on the right. Here we introduced the partial width Tg^ eq i = X' 2 Mg/167T 
for squarks undergoing $ p decays. It is seen that, in general, gauge decays contribute 
strongly at low Xi masses and small Yukawa couplings. 

The case A' 131 7^ (or A' 132 7^ 0) is of special interest || since it allows for direct 
production of the stop via e + d — > i (e + s — > t). The stop is particular in the sense that 
a "light" stop mass eigenstate (ij.) could (depending upon the mass parameters for the 
chiral states and on the free parameters of the model) exist much lighter than the top 
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Channel 



X? 
nature 



Decay processes 



Signature 



S5 



1,Z 



j,Z,H 



H 



Q Xi 
x' 

d xt 



vq'q' 



A' 



vq'q' 



q q 



d Xx 
x' 

d xt 



Missing Pt 
+ multiple jets 



vud 



•—^qq' 



S6 



H 



ul 



d xt 



W + X °i 
1+ v 



High Pt e + or /i + 
+ missing Pt + 1 jet 



S7 



1,Z 



d xt 



W + X ] 

^ e ± q'q" 
■ /+ v 



High P T e ± 
+ high Pt e + or fi^ 
+ missing Pt 
+ multiple jets 



SS 



d Xt 



W + X °i 

A' 



uq'q' 



/+ v 



High Pt e + or /i H 
+ missing Pt 
+ multiple jets 



Table 3: Squark decay channels in Ijt p SUSY classified per distinguishable event topologies 
(second part). As in table |2[ the list of processes given here is not exhaustive, e.g. the gauge 
decays xt ~^ Xi^ +l/ an d xt ~^ XiQQ r ma Y a ^ so proceed via a virtual sfermion. 



quark itself and lighter than other squarks. This applies only for the stop since the off- 
diagonal terms which appear in the mass matrix associated to the superpartners of chiral 
fermions are proportional to the partner fermion mass. Such a stop ti mass eigenstate 
is considered in this paper and its search is furthermore extended towards low mass by 
considering pair production via 7-gluon fusion as illustrated in Fig. ^. For the study of 
this process, we assume that the t\ is lighter than the lightest chargino. Hence the ti will 
decay dominantly into a positron and a quark since, by assumption, the decays into tx? 
and bxt are forbidden and the one-loop decay into cx\ is negligible even for small values 
of the coupling of the t\ to a positron-quark pair ||. 
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Figure 2: Stop pair production via 
7-gluon fusion at HERA, followed by 
Ijtp decay of the t\. 




3 The HI detector 

A detailed description of the HI detector can be found in ||. Here we describe only the 
components relevant for the present analysis in which the event final state involves either 
an e + (or e~) with high transverse energy or a large amount of hadronic transverse energy 
flow. 

The e + (or e~) energy and angle are measured in a finely segmented liquid argon 
(LAr) sampling calorimeter [ID| covering the polar angleQ range 4° < 9 < f53° and 



all azimuthal angles. It consists of a lead/argon electromagnetic section followed by 
a stainless-steel/argon hadronic section. Electromagnetic energies are measured with 
a resolution of a(E)/E ~ 12 %/y/E © 1% and hadronic energies with a(E)/E ~ 50 
%/VE © 2% after software energy weighting |TTJ. The absolute scales are known to 2% 



and 5% for electromagnetic and hadronic energies respectively. The angular resolution 
on the scattered electron measured from the electromagnetic shower in the calorimeter is 
< 4 mrad. A lead/scintillator electromagnetic backward calorimeter extends the coverage 
at larger angles (155° < 9 < 176°). 

Located inside the calorimeters is the tracking system used here to determine the 
interaction vertex and the charge of the final state lepton. The main components of this 
system are central drift and proportional chambers (25° < 9 < 155°), a forward track 
detector (7° < 9 < 25°) and backward proportional chambers (155° < 9 < 175°). The 
tracking chambers and calorimeters are surrounded by a superconducting solenoid coil 
providing a uniform field of 1.15 T within the tracking volume. The instrumented iron 
return yoke surrounding this coil is used to measure leakage of hadronic showers and to 
recognize muons. The luminosity is determined from the rate of the Bethe-Heitler process 
ep — > epy measured in a luminosity monitor. 



^^The incoming proton moves by definition in the forward (z > 0) direction with 6 = 0° polar angle. 
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4 Analysis 



4.1 Single production of squarks 

For the search for resonant production of squarks, the event selection basically relies on 
the final state lepton finding and on global energy-momentum conservation cuts. It is 
optimized separately for each of the event topologies (see tables || and ||) SI to S8 by 
relying on Monte Carlo simulation. 

The simulation of the leptoquark-like signatures (SI and S2) relies on the event 
generator LEGO |]r2[. For squarks undergoing gauge decays followed by a Xi or Xi Up 
decay into a high Pt e 1 * 1 and multiple jets, i.e. processes belonging to topologies S3 
and S4, the generator SUSSEX [JT^] based on the cross-sections given in 0] is used. 
Both generators also simulate initial state bremsstrahlung in the collinear approximation, 
initial and final state parton showers and fragmentation fl3], [14|], and properly take into 



account the correction of the kinematics at the decay vertex for effects of the parton 
shower masses. The parton densities used [15], |1| are evaluated at the scale of the new 



particle mass, and this scale is also chosen for the maximum virtuality of parton showers. 
For these channels, a complete simulation of the HI detector response is performed. The 
event topologies S5 to S8 (as well as some of the processes in S3 or S4 which proceed 
through the exchange of a virtual W or virtual sfermion), were studied at four- vector 
level || taking into account matrix element calculations [0 and multiparticle phase 
space. For these channels, realistic efficiencies are then obtained by smearing the particle 
four-vectors according to measured resolutions, detector effects and acceptances. The 
efficiencies thus obtained were cross-checked and found to agree typically within 5% with 
a complete simulation based on SUSSEX for those S5 processes where the xt undergoes 
a $p violating decay. 

A complete Monte Carlo simulation of the HI detector response is performed for each 
possible background source. For the DIS NC or CC background estimates we make use 



of either the DJANGO pf or the LEPTO f2(| event generator. DJANGO includes 



first order radiative corrections and simulation of real bremsstrahlung photons based on 



HERACLES as well as QCD dipole parton showers based on ARIADNE [[22]. LEPTO 
includes the lowest order electroweak scattering process with QCD corrections to first 
order in a s , complemented by leading-log parton showers and string fragmentation f23| . 
Both generators agree in channels where one expects a single hard jet, i.e. SI, S2 and 
S6. The LEPTO event generator is used in the multijet channels S3, S4, S5, S7 and 

S8. 

The parton densities in the proton used for DIS throughout are taken from the 
MRS H [[UJ parametrization which is close to F 2 structure function measurements at 
HERA (see ||24}|). For the direct and resolved photoproduction of light and heavy flavours, 



the PYTHIA MC event generator [|13| is used which includes QCD corrections to first or- 
der in a s , leading- log parton showers and string fragmentation [14[]- The GRV LO (GRV-G 
LO) parton densities |25| in the proton (photon) are used at low Q 2 . 



The event selection for real data starts with the rejection of non-colliding background. 
This selection step is common to all channels and requires: 
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1. a primary interaction vertex in the range | z — z |< 35 cm with z = 3.4 cm; 

2. that the event survives a set of halo and cosmic muon filters; for channel S2 these 
are complemented by visual scan; 

3. that the event be properly in time relative to interacting bunch crossings. 

Cut (1) mainly suppresses beam-wall, beam-residual gas and, with (2) and (3), back- 
ground from cosmic rays and halo muons. We moreover impose that the events be 
accepted by LAr calorimetry triggers the events of SI, S3 and S4 must satisfy 

"electron" or "transverse energy" trigger requirements; events of S2 and S5 must fulfill 
"missing transverse energy" requirements; events of S6, S7 and S8 must satisfy either 
"electron" or "missing transverse energy" requirements. 

The selection cuts and data reduction specific to each of the event topologies for the 
ep-induced background is presented below. In each case the number of event candidates 
observed are compared to SM expectations. The systematic (syst.) errors quoted on the 
mean expected background in each case take into account uncertainties on the absolute 
electromagnetic and hadronic energy scales (see section ^|), on the integrated luminosity 
(1.5%) and the contribution due to finite Monte Carlo statistics. Estimates of SUSY 
signal detection efficiencies are also given in each channel. 



Event topology SI: 

For the event topology SI, i.e. events characterized by the DIS NC-like signatures ||, 
it is necessary to reject contaminating background from other physical processes. We 
require: 



1. an isolated 'e + ' cluster |2§ with £ T>e = E e sin 6 e > 7 GeV and 10° <0 e < 145°; (here 
'e + ' includes all candidates except those having an associated track with explicitly 
measured negative charge); the isolation requires that less than 10% additional en- 



ergy be found within a pseudorapidity-azimuth cone of opening J (Arj e ) 2 + (A0 e ) 2 < 
0.25 centered on the e ± candidate; 

2. that if two 'e + ' cluster candidates are found, they must not be balanced in Ex, e and 
in azimuth, i.e. E\jE\ e > 1.25 and | A0 12 — 180° |> 2°, and the candidate with 
highest E Te must be at smallest rapidity; 



3. a total missing transverse momentum Px^miss ~ y(J2El,) 2 + (j2E^j < 15 GeV 
summed over all energy depositions i in the calorimeters, with E l x = E l sin Q % cos 0* 



and E % y = E % sin 9 l sin . 



4. a minimal "longitudinal momentum" loss in the direction of the incident positron 
—8 < 2E® — — P z ) < 12 GeV, where E® is the incident positron beam energy 

5. a y e , measured from the final state 'e + ', satisfying y e < 0.95. 
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Cuts (1) and (3) eliminate DIS CC events. Cut (2) suppresses QED Compton events. Cut 
(4) provides a powerful rejection of photoproduction contamination and also suppresses 
DIS NC-like events with a very hard 7 emitted from the initial state positron. Cut (5) 
further suppresses photoproduction with a "fake" e which tends to cluster at largest y e 
for largest M e , where y e is the standard DIS Lorentz invariant and M e the "squark mass" 
reconstructed from the final state 'e + ' energy E e and angle 9 e as : 

,2 E T,e = , _ E e - E e cos 6 e 

e ~l-2/ e ' Ve ~ 2E° 

Q 2 is the standard momentum transfer squared of DIS. In addition to the above require- 
ments, we apply a M e dependent y e cut which is designed |], |3[] via Monte Carlo studies 
to optimize the signal significance for scalar leptoquark searches, given the expected back- 
ground. This y e cut varies from y e > 0.5 at 45 GeV to y e > 0.35 at 150 GeV and down to 
y e > 0.05 at 275 GeV. 

For these NC-like (leptoquark like) signatures, 362 events satisfy the selection require- 
ments and the y e cut in the mass range M e > 25 GeV. This observed number of events 
is in good agreement with the mean expected DIS NC background of 335 ± 36 (syst.) 
events. The measured mass spectrum is compared to the DIS NC expectation in Fig. 0. 
For M e > 45 GeV, we are left with 91 events while 84 ± 10.2 (syst.) events are expected 





M e (GeV) M h (GeV) 



Figure 3: Mass spectra for (a) e + q (b) v + q final states for data (closed points) and 
DIS Monte Carlo (shaded histograms). The superimposed dashed histograms show typical JjL p 
SUSY signals near the sensitivity limit (see table [| in section ||) for (a) Mq = 150 GeV and 
M x o = 20 GeV and for (b) M q = 75 GeV and M x? = 20 GeV. 

from the SM. For M e > 100 GeV, 13 events are observed in good agreement with the 
mean SM expectation of 12.4 ± 2.6 (syst.). 
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In this channel, the $ p SUSY signal detection efficiency is found to be weakly depen- 
dent on Mg and ranges from 43% at 45 GeV to 68% at 150 GeV in the middle of the mass 
range considered here. 



Event topology S2: 

The event topology S2 is characterized by DIS CC-like signatures for which we require: 
1. no e ± cluster satisfying the above SI requirements; 
2- PT,miss > 25 GeV; 

3. the total transverse energy Et ~ 2 | Pt | calculated from energy depositions in the 
calorimeter should match the total missing transverse momentum Premiss such that 
(E T - P T ,miss)/E T < 0.5. 

Cuts (1) to (3) eliminate photoproduction and DIS NC background. 

In total, 40 CC-like events satisfy all above requirements in the relevant mass and y 
range at > 45 GeV and y^ < 0.95 where and yt are reconstructed by summing 
over all measured final state hadronic energy: 



Qh -o2 ^T.mias J2(P ~ P; 
> Qh - --> > Vh " 



\ Vh 1 - Vh 2E" 



As is seen in Fig. Qd, this is in good agreement with the DIS CC expectation of 33.4 ±3.6 
(syst.) events. 

The ty p SUSY signal detection efficiency in this channel rises from ~ 15% at 45 GeV 
to reach a plateau at ~ 80% above 100 GeV. 



Event topology S3: 

For a gauge decay of the squarks leading to a "right" sign final state lepton (i.e. e + q — > 
q — >• Xi + <?' ~~ > e + q"q"q'), we impose the following stringent requirements in complement 
to SI cuts : 

1. the 'e + ' must give y e > 0.4; 

2. an imbalance between the total hadronic E Tjh and Px t h such that (E Tjh —PT,h) / Er,h > 
0.25; 

3. at least one reconstructed jet with Prjet > 7 GeV; 

4. the azimuthal opening angle A<pj t between the jet of highest Pt (i.e. generally the 
current jet in a DIS NC process) and the axis defined by the total hadronic transverse 
momentum P^/i = {J2 E Xt h,J2 E y ^) be larger than A(f>j t > 2/5 x (50 — Et^) with 
Et,h in GeV and Acj)j t in degrees; 
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5. the squark invariant mass calculated from all final state particles (M<fe c ) excluding 
the proton fragments [1] must deviate from M e by more than 10%. 

Cut (1) strongly suppresses DIS NC background. The y e distribution calculated from 
the 'e + ' in SUSY events appears strongly shifted towards large y e since the squark decays 
uniformly in its center-of-mass frame and, further, since the 'e + ' takes away only a fraction 
of the Xi momentum. Cut (2) exploits the fact that the hadronic energy of the event is not 
concentrated within one single jet because of the decay products of the Xv Fig. [|a shows 
how this value discriminates the DIS NC background from the signal. The jet finding for 
cuts (3) and (4) relies (here and throughout the paper) on a simple cone algorithm in the 
laboratory reference frame with a fixed pseudorapidity- azimuth opening radius of 1 unit. 
Cut (4) further suppresses lowest order DIS NC events by imposing sufficient hadronic 
activity far enough in azimuth from the "current" jet. Cut (5) ensures that the events 
accepted here (gauge decay modes) are not simultaneously accepted as SI candidates ($ p 
decay modes). 




Figure 4: (a) Balance between the total scalar hadronic energy Eth an d the vector sum Pt,h 
for DIS NC before and after the y e > 0.4 cut and for a simulation of a 150 GeV squark decaying 
into e + q — ► e + q"q"q'; (b) measured mass spectrum (uncorrected) for the squark candidates 
in the e + q — ► e + q"q"q' channel for data (dots) and DIS NC Monte Carlo (histogram), the 
superimposed dashed histogram in (b) show a typical Il p SUSY signal near the sensitivity limit 
(see table | in section §) for Mg = 150 GeV and M x o = 40 GeV. 

We find 405 candidates satisfying the previous cuts for masses M^ ec > 25 GeV, which 
is to be compared with the mean SM background of 363±39 events expected from DIS NC. 
The measured mass spectrum is compared to Monte Carlo expectations in Fig. [|b. Above 
45 GeV we observe 220 events in the data, while 154 ± 17 (syst.) events are expected 
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from DIS NC. This represents an excess of 2.9 standard deviations in the Gaussian limit 
approximation (combining statistical and systematic errors in quadrature). The slight 
excess of events is seen to be mostly concentrated at low masses and in particular in the 
mass range from 40 to 85 GeV. Nevertheless, it should be recalled here that our DIS NC 
Monte Carlo for multijet channels (LEPTO) does not include full QED corrections which 
could lead to a migration of events with true y e below cut (1) towards larger apparent 
y e . A good agreement is observed for M dec > 100 GeV where we find 14 events while the 
mean DIS NC expectation is 13.2 ± 2.6 (syst.) events. For M dec > 140 GeV, 5 events are 
found which agrees well with the expectation of 1.9 ± 0.9 (syst.). 

In this channel, the E\ v SUSY signal detection efficiencies (which experimentally sums 
that of both the right sign and the unsigned events), depend mainly on the Xv For 
M x o = 20 GeV it rises from ~ 20% for M q = 45 GeV to a plateau of ~ 33% for M q > 
75 GeV. For M x o = 80 GeV it reaches ~ 60% for M q > 100 GeV. 



Event topology S4: 

For a gauge decay of the squarks leading to a "wrong" sign final state lepton (i.e. e + q — > 
q — > e~q"q"q') we perform a determination of the lepton charge using the tracking chamber 
information. Hence, we impose in addition to the above S3 criteria : 

1. the e~ LAr cluster must be geometrically linked to a negatively charged track and 
the cluster energy must match the track momentum within | (E — P)/(E+P) |< 0.5; 

2. the track must be made of at least 40 digitisations in the central tracking chamber; 

3. the error in the curvature k must fulfil | k/8k |> 1. 

These cuts ensure a good quality of the track reconstruction and track-cluster matching 
at the expense of a reduced efficiency (partly due to occasional inoperation of either the 
inner or the outer central drift chambers) for accepting the e~ track of about 70% in the 
angular range 9 > 35° well covered by the central tracking chambers. 

We observe no e~ (wrong sign) events among the 405 candidates satisfying the kine- 
matical requirements for squark gauge decays. 

In this channel, the charge track requirements imply an additional efficiency loss com- 
pared to S3 for the I/t p SUSY signal which is negligible at 45 GeV but which increases to 
10% at 150 GeV and 20% at 250 GeV. 



Event topology S5: 

The event topology S5 is characterized by large missing transverse momentum and mul- 
tiple jets. We require: 

1. no e ± cluster satisfying the above SI requirements; 

2. P T ,miss > 15 GeV; 
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3. {E T>h - P T , h )/E Tth > 0.25; 

4. P TA > 50 x (1 - (E T>h - P T>h )/E T>h ) with P TA in GeV; 

5. at least one reconstructed jet with Pr,jet > 7 GeV; the jet of highest Pt should 
satisfy A(fij t > (4/7) x (100 — Er,h) with A(f)j t in degrees and Er,h in GeV. 

Cut (3) exploits the fact that more than one jet is expected in such events. Cut (4) removes 
the DIS CC background, which is mainly concentrated at low values of (E T)h — P T h )/E T>h . 
Cut (5) removes photoproduction events for which one of the two back-to-back jets is badly 
measured, so that the is expected to be small. We are left with 9 events in the data 
sample compared to an expectation of 3.9 ± 4 (syst.) events from jp photoproduction 
background and a negligible DIS CC background. Here a sizeable contribution to the 
systematic error originates from the dependence of the LAr trigger efficiency on this S5 
multijet topology. 

It is shown in Fig. ||how these 9 remaining events compare to the SM photoproduction 
and DIS CC expectations (respectively for 1 x Cdata and 10 x Cdata)- No LAr trigger 
efficiency losses are folded in the Monte Carlo sample of Figs. |5|b,c, and d. It can be seen 
also that cut (5) still ensures a good efficiency for a possible SUSY signal. 

The efficiency for $ p SUSY events in S5 rises with increasing Mq up to a plateau for 
Mq > 150 GeV. It also rises with increasing M x o mainly because of the PT,miss selection 
cut imposed. For processes where the final state Xi ( or Xi) at the squark decay vertex 
directly undergoes a$ p decay, the efficiency at Mq > 150 GeV is ~ 26% for M % o = 20 GeV 
and ~ 52% for M x o = 80 GeV. When the Xi is H-\ike (stable) the mass difference between 
the Xi an d the xt * s relatively smaller for heavier x± and this hampers the signal detection. 
Hence in processes where the Xi escapes detection, the efficiency which at Mq > 150 GeV 
reaches ~ 50% for M x o = 20 GeV is down to ~ 30% for M x o = 80 GeV. Finally, processes 
where the xt undergoes a cascade decay (e.g. xf ~^ W + x± ', Xi ~^ e + q'q" ; W + qq'), 
suffers from an efficiency loss due to the PT,miss cut. In such case the efficiency rises from 
~ 6% at M q ~ 100 GeV to ~ 32% at Mq > 175 GeV. 



Event topology S6: 

The event topology S6 is characterized by the presence of a lepton (e + or at large 
transverse energy Etj, a large missing transverse momentum and a single jet. 
To search for cases where the final state lepton is a positron, we require: 

1. an isolated 'e + ' with Et^ > 7 GeV; 

2. P Timiss > 15 GeV; 

3. (E Tth -P Tth )/E Tjh <0.B; 

4. 0.4 <y e < 0.95; 

5. at least one reconstructed jet with Prjet > 7 GeV; 
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Figure 5: Correlation between the azimuthal opening angle Afyt and the total hadronic 
transverse energy £r,ft for (a) data, (b) DIS CC and background MC and (c), (d) Ijt p SUSY 
signal in topology S5 for two example cases. The events above the cut (solid line) are accepted. 

6. at least 1 charged track with Pt rac k > 5 GeV originating from the primary vertex 
and linked to the electron cluster. 

Cut (1) suppresses DIS CC background while cut (2) suppresses photoproduction back- 
ground. The other cuts are designed to optimize the specific S6 signal significance relative 
to tails of background distributions. We are left with 2 event candidates in the data while 
3.8 ± 1.3 (syst.) are expected from DIS NC background. 
To search for cases where the final state lepton is a muon, we require: 

1. no isolated e ± with E T)( , > 7 GeV; 

2- PT,miss > 25 GeV calculated from the energy depositions in the calorimeters; 
3. (E TA - P T , h )/E Tjh < 0.5; 
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4. at least one reconstructed jet with Prjet > 7 GeV; 

5. at least 1 charged track with Ptrack > 10 GeV linked to the primary vertex and lying 
outside of a A0 = 60° cone centered on the direction of the jet at highest Prjet] 

6. that the charged track in cut (5) be "penetrating" in the sense that there be < 5 GeV 
of total energy measured in the LAr calorimeter within a A0 =15° cone centered 
on the track direction; moreover the track should not point to a cluster localized 
around the azimuthal cracks in the LAr calorimeter. 

We are left with 1 event candidate in the data of which an event display is shown Fig. ^|. A 
remarkable "/i + jet" signature is seen (Fig. |6|a) with a positively charged isolated track. 
A detailed analysis of this event [^] reveals that the isolated track has a transverse 




c) 

Figure 6: Display of the e + p — > fi + jet candidate from HI run number 84295 showing (a) 
an R-z view, (b) a H-4> view and (c) the calorimetric transverse energy flow. 

momentum of 23 ± 2A+1 GeV and appears in azimuth at A0 M = 183 ± 1° from an 
hadronic system (Fig. ^o) which itself has in total Px,h = 42.1 ± 4.2 GeV. This hadronic 
system is built from two main "clusters" (Fig. |]c) which are found to be merged into a 
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single jet with our cone algorithm within a radius of 1 unit in the pseudorapidity-azimuth 
plane. This leaves overall a total missing transverse momentum of | Pt,k + Pt,h |— 
18.7 ± 4. 8^7 GeV and a "longitudinal momentum" loss in the incident positron direction 
of 2E° e - {£ (E h - P g>h ) + - P^)} = 35.8 ± 1.61^ GeV. 

With our selection cuts (1) — > (6), we expect on average from DIS CC a background 
of 0.04 ± 0.03 (syst.) event corresponding to a probability of about 4% to observe here 
one or more such event. The actual dominating background is expected |27|] to come 
from associated W + production in NC and CC processes followed by a leptonic decay of 
the W + , e + p — ■> e + W + X — > e + fi + uX, where the final state e + is lost in the beam pipe. 
With our selection cuts, this process studied in a Monte Carlo calculation using the SM 
cross-section of |28j gives a 15% probability for such a background event to be accepted. 

The efficiency for J^ v SUSY events in S6 rises with increasing Mg up to a plateau 
for Mg > 150 GeV, and decreases with increasing M x o for reasons already explained in 

topology S5 for P-like xl- For Mg > 150 GeV it is of ~ 50% for M x o = 20 GeV, and 
down to ~ 35% for M Y o = 80 GeV. 

It is interesting to note [gij that an E\ v SUSY signal is consistent with the properties 
of the observed event candidate. Imposing more restrictive requirements such as: 2E® — 
{E(E h -P zA ) + (E^-P^)} > 12 GeV, A0 M > 140°, P T , h > 40 GeV and P T , M > 
10 GeV, we find that more than 30% of the l?t p SUSY events satisfying the selection 
cuts (1) to (6) in the muon channel also verify these additional requirements provided 
that the ul be in the mass range 100 < M^ L < 200 GeV and the xf be in the range 
20 < M v+ < 30 GeV for M Sr ~ 100 GeV or 20 < Af + < 110 GeV for M Sr ~ 200 GeV. 

With the above stringent requirement that Pt,h > 40 GeV applied on the background, 
which implies in associated W + production a very stiff recoiling hadronic system, one 
is left with a ~ 3% probability for such interpretation of the event and a negligible 
contribution from misidentified DIS CC events. 



Event topology S7: 

The event topology S7 is characterized by the presence of an e ± at large E T , accompanied 
by another lepton l + , large missing transverse momentum and multiple jets. We require: 

1. an isolated e ± with Et^ > 7 GeV and giving 0.4 < y e < 0.95; 

2. P Timiss > 15 GeV; 

3. (Eta. - P T ,h)/E T<h > 0.25; 

4. at least one reconstructed jet with Prjet > 7 GeV. 

In the case where Z + is a positron we require in addition a second isolated e with E? > 
5 GeV. We observe no candidate in the data while 0.4 ± 0.3 event is expected on average 
from DIS NC and none from DIS CC. In the case where Z + is a muon we impose a more 
stringent cut on PT,miss > 25 GeV. We observe no candidate in the data, while 0.4 ± 0.3 
event is expected on average from DIS NC and none from DIS CC. 
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The efficiency for $ p SUSY events in S7 depends weakly on Mg and rises with M x o. 
It is of ~ 20% for M x o = 20 GeV and rises in the case where l + is a positron to ~ 60% 
for M x o = 80 GeV. In the case where l + is a muon and for heavy Xi, the apparent Premiss 
tends to be reduced by the quasi-collinearity of the \i and neutrino in the final state and 
the efficiency only reaches ~ 30% for M o = 80 GeV. 



Event topology S8: 

The event topology S8 is characterized by the presence of one charged lepton (e + or 
at large Et, a large missing transverse momentum and multiple jets. 
For the case where the lepton is a positron, we require: 

1. an isolated 'e + ' with E T , e > 7 GeV and giving 0.4 < y e < 0.95; 

2. P Timiss > 15 GeV; 

3. (E Tih - P Tth )/E TA > 0.25; 

4. at least one reconstructed jet with Prjet > 7 GeV. 

We are left with 3 candidates in the data while 2.3 ± 1.0 are expected from DIS NC. The 
background from DIS CC and photoproduction is here negligible. 
For the case where the lepton is a muon, we require: 

1. no isolated with Et^ > 7 GeV; 

2- Premiss > 25 GeV; 

3. (E T>h - P T , h )/E Tth > 0.25; 

4. at least one reconstructed jet with Prjet > 7 GeV; 



5. at least 1 "penetrating" and isolated charged track with P tra ck > 10 GeV, such 
tracks were defined in cuts (5) and (6) of the S6 selection. 

We are left in the data with the same /i + X event candidate as was found in S6 event 
topology while 0.04 ± 0.03 event is expected from DIS CC background. 

The efficiency for fi p SUSY events in S8 has similar Mg and M x o dependence as in 
S7. It is of ~ 20 — 30% for M x o = 20 GeV and rises in the case where Z + is a positron 
to ~ 65% for M x o = 80 GeV. In the case where l + is a muon, it remains at ~ 25% for 
M v o = 80 GeV. 

While the efficiency is rather high for the ty p SUSY signal to satisfy the basic cuts 
(1) to (5), it is here (contrary to the S6 channel) difficult to meet the more restrictive 
conditions |29| that could be imposed on the observed fi + X event candidate. In particular, 



the overlap of hadronic Pt flow initiated by the three quark jets into a single observed 
jet and with A^^^ > 140° is a unlikely configuration. It is moreover difficult to satisfy 
simultaneously the stringent cut 2£° - (E h - P Zjh ) + (E fl - P Z4l )} > 12 GeV. We find 
that less than 10% of the ljt p SUSY events satisfying the selection cuts (1) to (5) also 
verify these additional requirements for be in the mass range 100 < M^ L < 200 GeV. 
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4.2 Pair production of stop squarks 

For squarks which are pair produced in 7-gluon fusion processes (Fig. the scattered 
electron is generally lost in the beam pipe. The optimization of the event selection relies 
on Monte Carlo simulation. The simulation of stop pair production in 7-gluon fusion |J is 
based on the a stop cross-section calculated in the Weizsacker- Williams approximation. The 
resolved photon contribution |34| as well as the contribution from a Zq boson exchange 



are neglected. The GRV LO gluon density in the proton [25| is used. The background 
simulation is based on the event generator DJANGO H for DIS NC and PYTHIA @ 
for photoproduction processes. We impose: 

1. P T , m iss < 15 GeV; 

2. two candidates i satisfying E l T e > 5 GeV within 10° < 9\ < 145°, the must 

be isolated within pseudorapidity-azimuth cones of opening yj (Ar/*) 2 + (A</>*) 2 < Ri 
where Ri=i = 0.5 and Ri=2 = 0.25; 

3. the two candidates must be acollinear in the transverse plane, A0 12 < 140°; 

4. there must be missing "longitudinal momenta" such that 2E® — — P z ) > 
12 GeV; 

5. at least two jets must be found by the jet cone algorithm each with Exjet > 5 GeV. 

Cut (1) suppresses the contamination of bb photoproduction where the heavy b quarks can 
undergo semi-leptonic decays. The transverse energy requirement in cut (2) is high enough 
to eliminate contamination from J/ip photoproduction. Asking for two isolated e ± strongly 
suppresses the main DIS NC and photoproduction background, leaving only events where 
either a photon or a hadronic jet is misidentified as an e ± . Such a misidentified jet in DIS 
NC events is most probably found to be collinear with the true electron in the transverse 
plane. This can be seen in Fig. [f]a where the azimuthal balance A0i t 2 between the two 
candidates found in the data is compared with SM Monte Carlo simulation. For this 
comparison with DIS NC, an additional cut of J2 (E — P z ) > 43 GeV is imposed. The 
observed acollinearity in between the two e ± candidates is seen to be well described by 
the Monte Carlo. 

The remaining background is suppressed with cut (3) and we are left with 26.1 ± 3.9 
expected events from DIS NC while the photoproduction contamination is negligible. The 
estimated DIS NC background agrees well with the 28 observed events. This background is 
further suppressed by cut (4). The effect of this cut is seen in Fig. 0b showing Monte Carlo 
expectations for the correlation between ~ Pz) and the azimuthal balance A0 12 for 
DIS NC (plotted for 3 x Cdata) and stop pair production (arbitrary normalization). Two 
events are rejected by cut (5). 

With cuts (1) to (5), we are left with no observed event candidate while 1.0 ± 0.8 
event is expected from misidentified DIS NC background. The detection efficiencies for 
the signal of ti pair production after all selection cuts are found to rise from about 3% at 
9 GeV to 32% at 24 GeV. 
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Figure 7: (a) Azimuthal balance A0i,2 of the two e 1 * 1 candidates for data (closed dots) and 
DIS Neutral Current Monte Carlo (histogram) when {E ~ Pz) > 43 GeV. (b) Comparison 
of Monte Carlo expectations for the missing longitudinal momentum (E — P z ) versus the 
balance in azimuth of the two candidates Ac/>i,2 for DIS Neutral Current (closed dots) and 
stop pair production in 7-gluon fusion (open squares) for Mf = 20 GeV. 



5 Results 

In the absence of significant deviation from the SM expectations, we now derive exclusion 
limits for the Yukawa couplings \[j k as a function of mass, combining all contributing 
channels and making use of the number of observed events, of expected background events, 
and the signal detection efficiencies for each contributing channel. 



Single production of squarks: 

For the event topologies SI to S3, the detection efficiencies are folded with a mass bin 
of variable width which slides over the accessible Mg range. The bin width is optimized 
taking into account the mean expected background and the expected q mass resolution 
and contains about 68% of the signal at a given mass, e.g. at 150 GeV we typically have a 
full bin width of AM e ~ 25 GeV in SI, of AM h ~ 35 GeV in S2, and of AM dec ~ 40 GeV 
in S3. For event topologies S4 to S8 where the number of observed events and expected 
background is always < 10, the signal is integrated above the selection cuts (i.e. without 
explicit restriction on the reconstructed mass). More details on the methodology for the 
limits derivation and on the procedure for folding the channel per channel statistical and 
systematic errors are given in p]. 

The detection efficiencies have to be folded with the branching fractions in each of 
the possible event topologies, properly taking into account the relative production cross- 
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section of the various squark flavours (i.e. da and ul for A' m 7^ 0). For small M x o/Mq 
and/or small X[j k values where gauge decays of the squark are expected to dominate, 
the dependence upon the values of the free parameters of the MSSM has to be fully 
considered for the coupling constants at the q — > q + Xi or Q ~~ * <?' + Xi vertex as well 
as the decay branchings of the Xi an d xf- We will discuss our results in terms of the 
usual parameters: the ratio tan (3 of the two Higgs field vacuum expectation values, the 
higgsino mixing parameter fi, the mass parameter M2 for the SU(2) gauginos. 

The domains of the (M2,/!) plane where decays into a specific neutralino or chargino 
dominate for the Ul are shown in Fig.^a for a squark mass in the middle of the accessible 
Mg range. The decay into the lightest chargino xt is seen to dominate as soon as kinemat- 
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Figure 8: (a) Regions in the M2, [i plane where a 150 GeV ul squark decays dominantly into 
one of the %i or xf states for tan/3 = 1; (b) Main component of the lightest neutralino in the 
plane (M2,//), for tan/3 = 1. In (a) and (b) the region below the full curve corresponds to a 
domain in which the branching ratio of the Xi i n to e^qq' is greater than 80%. Along the dotted 
curve, the branching into e^qq' is of about 20% for a Xi °f 20 GeV. 



ically allowed. This is contrary to the cIr which mainly decays into d + Xx (except when 
the xi is H- dominant). For a if-dominant LSP the coupling constant is small enough to 
suppress strongly squark gauge decays into Xi- 

Regions for various kinds of Xi m the parameter space are shown in Fig. [8]b. Similar 
plots of the admixture of the weak eigenstates 7, Z, H® and H® in the neutralino LSP are 
given in [32]. In the region where the Xi is 7-like, the branching ratio of the Xi i n l° e^qq' 
is greater than w 60% for tan/3 = 1. When the Xi is a pure 7 (M 2 = 0), this branching 
saturates at about 88%. In the region where the x% is Z-\ike, for instance along the dotted 
line for a Xi °f 20 GeV, the branching ratio into e^qq' is about 20% 
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Table 4: Number of observed events, of expected background events and branching ratios 
B (in %) corresponding to channels SI to S8 for some values of Mq and M x o, and presented 

for 7-like, Z-like and iJ-like Xi- The efficiencies e are given in each channel for the 7-like case 
except* for S6 where they are given for the ii-like case. Also given is the J2 £ & summed over 
all channels for each Xi type. 
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The relative contribution of each of the channels SI to S8 is given in Table [| for a 
few representative cases. For the chosen examples, the mass of the xf is about twice that 
of the Xi expect for heavy (e.g. 80 GeV) H-like xi f° r which M x +/M x o ~ 1.3. In the low 
squark mass region, the decay modes S3 and S4 dominate since the decay of the squark 
into xt is largely suppressed by phase space and since the Yukawa couplings probed are 
small. For medium masses, when the q is heavier than the xi: it decays dominantly via 
S3 and S5. The channel S5 contributes mainly when the xi is -Z-like, since in that case 
the branching ratio for x± v + 2jets is generally above 70%. At very large q masses, we 
are only sensitive to high values of the coupling, so that the $ p decay SI into e + + jet 
dominates whereas ftp decays S2 are strongly suppressed by the parton density. 

From the total branching fraction (i.e. summed over all above channels), it is inferred 
that the contribution of the decays into heavier x°(i > 1) and xf (j > 1) is generally 
small. Hence, in order to simplify the derivation of limits we have assumed conservatively 
that the q decay into these heavy states are allowed but measured with vanishingly small 
efficiencies. Folded in the derivation of limits are systematic errors coming from the 
uncertainty on the luminosity measurement (1.5%), the finite Monte Carlo statistics, the 
absolute energy calibration which leads to an uncertainty on the background estimation 
of about 10%, the choice of the scale entering the structure function calculation (which 
leads to « 7% uncertainties in the cross-section) and the choice of the parton density 
parametrization. 

The exclusion limits obtained on the coupling A' m at 95% confidence level (CL) are 
shown in Fig. ^| as function of the q mass (ul and da) in the hypothesis that the Xi is 
a pure 7. The limits represent an improvement of about a factor two to three compared 
to our previously published results [0] at low squark masses. The gain from integrated 
luminosity is only partially cancelled by a less favourable quark density since for incident 
positron, the e + d — > ut production dominates over e + u — *> d^ whilst e~u — > d^ dominates 
over e~d — *> ul for incident electron data. The limits also improve at largest mass where 
the smaller coupling probed implies a narrow observable resonance width. 

For the smallest couplings accessible here, a squark lighter than m 230 GeV undergoes 
dominantly a gauge decay. Hence in contrast to earlier searches M , we are here sensitive 
to event topologies immediately distinguishable from those of leptoquarks. The existence 
of first generation squarks with fi p Yukawa coupling A' m is excluded for masses up to 
240 GeV (depending on the x± mass) at coupling strengths A' 1 2 11 /47r > a em ( up to 130 GeV 
for coupling strengths > 0.01 x a em ). 

In the more general case, where the xi is a mixture of gauginos and higgsinos, the 
exclusion limits on the coupling A' m at 95% CL are shown in Fig. [T^ as function of 
the q mass. Here, the limits are derived for a reference point in the MSSM parameter 
space chosen as /i = —160 GeV and M 2 = 60 GeV for a 7-like xi, A 4 — 150 GeV and 
M 2 = 150 GeV for a Z-like xi, and // = -44 GeV and M 2 = 140 GeV for a H-hke x°v 
These points in the parameter space lie outside the domain excluded from the invisible 
Zq width measurement at LEP p3| . The rejection limits for a 7-like Xi are found not to 
differ much from those obtained for pure 7. The limits are also seen not to depend too 
strongly on the nature of the LSP. The three curves merge together at highest squark 
masses, where the branching of the squark into x?<? becomes negligible relative to $ p 
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Figure 9: Exclusion upper limits at 95% CL for the coupling A' m as a function of the squark 
mass for various fixed photino masses and derived for tan/3 = 1. The regions above the curves 
are excluded. The limits combine I/t p and gauge decays of the cIr and u^. 

decay. The limits obtained are found moreover not to depend strongly on the parameter 
tan f3. Varying tan f3 from 1 to 40, we find that the limits only slightly degrade and mainly 
at very low q masses, by 30% at Mq = 45 GeV down to 2% at 200 GeV. 

From the analysis of the A' m case involving the cLr and ut squarks, limits can be 
deduced on the A' ljfc by folding in the proper parton densities. Such limits are given in 
Table ^| at Mq = 150 GeV. For Mq > 150 GeV, the exclusion limits for A' m and A' 121 in 
particular are found to coincide within 5%. 

Our rejection limits extend considerably beyond the only other collider limits of 
Mq > 100 GeV for A' m inferred in [[3(J from dilepton data of the Tevatron experiments. 
Moreover, this Tevatron limit was derived only under the restrictive assumption that the 
LSP is the 7 and that squarks other than the stop are degenerate in mass. If one assumes 
that one squark is substantially lighter than the others, this bound is much weaker || and 
only slightly above the mass reach of LEP 1. On the contrary, our A' m coupling limits 
only weakly depend on the mass degeneracy assumption since, as mentioned above, the 
u production strongly dominates. Hence, if the first generation u is significantly lighter 
than other squarks, we probe here (and similarly for the d in [1]]) a large portion of the 
mass-coupling domain unexplored by other experiments. 
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Figure 10: Exclusion upper limits at 95% CL for the coupling A' m as a function of squark 
mass for different natures of the x? f° r = 40 GeV and for tan (3 = 1 (region above the 
curve excluded). The full curve corresponds to cases where the LSP is 7-like, the dotted one to 
a //-dominant LSP and the dashed one to a Z-like LSP. 

There are no other direct limits published for \' x - k (j or k 7^ 1). This is particularly 
interesting since indirect limits for this coupling are also weaker H than those []35]| for 

^111 ■ 

Production of stop squarks: 

For the derivation of exclusion limits for the pair production of stop squarks, a 15% 
uncertainty on the cross-section (which varies from a s t op ~ 200 pb at 9 GeV to a stop = 
1 pb at 24 GeV) is due to the specific choice of gluon density. The uncertainty was 
determined by comparing with MRSD- [T5| and constitutes the main source of systematic 



error. 

The exclusion limit obtained on the stop mass at 95% CL is shown in Fig. [11]. A stop 
in the range 9 < < 24.4 GeV is excluded at 95% CL. This limit does not depend 
upon the value of \' 13k (as far as A' 13fc x cos 6^ > 10~ 4 , below which value the decay into 
c and Xi cannot be neglected anymore). The angle 9t is the mixing angle in the mass 
matrix of the stop (see for example |36| ) . 

If the coupling A' 131 dominates, the stop can be singly produced in reactions of the 
type e + + d — > t. Under our phenomenological assumptions (see section ^) and as long 
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Table 5: Exclusion upper limits 
at 95% CL on the couplings A' lj7c for 
Mg = 150 GeV and M %0 = 80 GeV. 
The quoted values for A' m are given 
for 7-dominant and for Z-dominant 
Xi- In other cases, the higgsino com- 
ponent of the Xi is assumed to be van- 
ishingly small. Moreover, the results 
for cases with j = 3 are only valid 
under the additional restriction that 
M x+ > M q . 
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as (A' 131 x cos 9 t ) > 10 -4 , the search of the t borrows from the analysis for ty p decays of 
first generation squarks. We find that masses below 138 GeV are excluded at 95% CL 
for coupling strength of (A' 131 x cos^t) 2 /^ > 0.01a em . This represents an increase of 
sensitivity of about an order of magnitude compared to our previous results [0]. As a 
comparison with other experiments, a t lighter than 38 GeV is excluded at 95% CL from 
LEP data [[37| for 9 t = from the width of the Zq. But for 9 t close to the value for which 
the Zq decouples from the stop (0.7 < 9 t < 1.4) there are no existing limit from LEP for 
Ijtp stops. 

The coupling limit in Fig. 11 is extended beyond Mi ~ M top . This portion of the 
exclusion limit curve is only valid for Mi < M top + M x o . 



6 Conclusions 

We have searched for squarks from i?-parity violating supersymmetry. The search was car- 
ried out for the first time at HERA in all possible decay processes allowed when wandering 
in the parameter space of the Minimal Supersymmetric Model. No significant evidence 
for the production of squarks was found and mass dependent limits on the couplings were 
derived. The existence of first generation squarks at masses up to 240 GeV are excluded 
at 95% confidence level for a strength of the Yukawa coupling A' m of A' 1 2 11 /47r = a em . The 
limits extend far beyond results obtained at other colliders where our excluded domain in 
the mass-coupling plane for masses > 100 GeV has never been explored. 

Scalar stop squarks were searched in pair and single production modes. The existence 
of light scalar stops with A' 13fc couplings to light fermions is excluded for masses 9 < mi < 
24.4 GeV at 95% confidence level. Stop squarks with A' 131 couplings are excluded below 
138 GeV at 95% confidence level for couplings (A' 131 x cos# 4 ) 2 /47r > 0.01a em . 
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Figure 11: Exclusion limits for the coupling A' 131 x cos9 t as function of the stop mass. The 
exclusion curve is obtained from the single stop production (region above the curve excluded). 
Beyond the dashed line, the limits are only valid for Mf < Mt op + M^o . The hatched domain is 
excluded by the pair production search and concerns A' 13fc with k = 1,2. 

Acknowledgements 

We wish to thank the HERA machine group as well as the HI engineers and technicians 
who constructed and maintained the detector for their outstanding efforts. We thank the 
funding agencies for their financial support. We wish to thank the DESY directorate for 
the support and hospitality extended to the non-DESY members of the collaboration. 

References 

[1] HI Collaboration, T. Ahmed et al., Z. Phys. C64 (1994) 545. 

[2] J. Butterworth and H. Dreiner, Nucl. Phys. B397 (1993) 3, and references therein. 

[3] HI Collaboration, S. Aid et al, Phys. Lett. B369 (1996) 173. 

[4] W. Buchmiiller, R. Riickl and D. Wyler, Phys. Lett. B191 (1987) 442. 



29 



[5] E. Perez and Y. Sirois, "SUSY Searches at HERA" , To be published in the Proceed- 
ings of the International Workshop on Supersymmetry and Unification of Fundamen- 
tal Interactions, Editions Frontieres (I. Antoniadis and H. Videau, Editors), (May 
1995, Palaiseau). 

[6] E. Perez, "Recherche de Particules en Supersymetrie Violant la R-parite dans HI a 
HERA", These de Doctorat, DAPNIA/SPP report (in French), to be published. 

[7] J.F. Gunion and H.E. Haber, Nucl. Phys. B272 (1986) 1. 

[8] T. Kon, T. Kobayashi and S. Kitamura, Phys. Lett. B333 (1994) 263; T. Kon et al., 
Z. Phys. C61 (1994) 239. 

[9] HI Collaboration, I. Abt et al, "The HI Detector at HERA", DESY preprint 93-103 
(July 1993) 194pp; idem DESY Internal Report Hl-96-01 (March 1996) 157pp. 

10] HI Calorimeter Group, B. Andrieu et al., Nucl. Instr. and Meth. A336 (1993) 460. 

11] HI Calorimeter Group, B. Andrieu et al, Nucl. Instr. and Meth. A350 (1994) 57; 
idem, Nucl. Instr. and Meth. A336 (1993) 499. 

12] LEGO 0.02 and SUSSEX 1.5; K. Rosenbauer, Doktor Dissertation, RWTH Aachen 
(in German), PITHA 95/16 (July 1995). 

13] PYTHIA 5.6; T. Sjostrand, Comp. Phys. Comm. 39 (1986) 347; 
T. Sjostrand and M. Bengtsson, Comp. Phys. Comm. 43 (1987) 367. 

1] JETSET 7.3; T. Sjostrand, CERN preprint TH-6488-92 (May 1992) 284pp. 

5] A.D. Martin, W.J. Stirling and R.G. Roberts, Phys. Rev. D47 (1993) 867. 

3] A.D. Martin, R.G. Roberts and W.J. Stirling, Durham Univ. preprint DTP-93-86 
and Rutherford Appleton Lab. preprint RAL-93-077 (October 1993) 16pp.; 
(PDFLIB |T7J nucleon structure function type 1, group 3, set 36) . 

7] H. Plothom-Besch, CERN-PPE Parton Density Functions program W5051. 

J] A. Bartl, H. Fraas and W. Majerotto, Z. Phys. C30 (1986) 441. 

3] DJANGO 2.1; G.A. Schuler and H. Spiesberger, Proceedings of the Workshop 
Physics at HERA, W. Buchmuller and G. Ingelman (Editors), (October 1991, DESY- 
Hamburg), vol. 3 p. 1419. 

[20] LEPTO 6.1; G. Ingelman, Proceedings of the Workshop Physics at HERA, 
W. Buchmuller and G. Ingelman (Editors), (October 1991, DESY-Hamburg), vol. 
3 p. 1366. 

[21] HERACLES 4.4; A. Kwiatkowski, H. Spiesberger and H.-J. Mohring, Corn- 
put. Phys. Commun. 69 (1992) 155. 

[22] ARIADNE 4.0; L. Ldnnblad, Comput. Phys. Commun. 71 (1992) 15. 

[23] JETSET 7.4; T. Sjostrand, Lund Univ. preprint LU-TP-95-20 (August 1995) 321pp; 
idem, CERN preprint TH-7112-93 (February 1994) 305pp. 

[24] HI Collaboration, I. Abt et al., Nucl. Phys. B407 (1993) 515, 

ZEUS Collaboration, M. Derrick et al, Phys. Lett. B316 (1993) 412. 



30 



[25] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D45 (1992) 3986; idem, Phys. Rev. D46 
(1992) 1973; (PDFLIB |HJ nucleon structure function type 1, group 5 set 4; photon 
structure function type 3 group 5 set 3). 

[26] HI Calorimeter Group, B. Andrieu et al, Nucl. Inst, and Meth. A344 (1994) 492. 

[27] HI Collaboration, T. Ahmed et al., "Observation of an e + p — ► fi + X Event with High 
Transverse Momenta at HERA" , DESY preprint 94-248 (December 1994) 9pp. 

[28] U. Baur, J.A.M. Vermaseren, D. Zeppenfeld, Nucl. Phys. B375 (1992) 3. 

[29] In the course of completing this paper, we became aware of an alternative ljt v SUSY 
interpretation of the observed e + p — > fi + X event candidate proposed by 
T. Kon, T. Kobayashi and S. Kitamurama, Seikei Univ. preprint ITP-SU-96/02 
(Tokyo, January 1996) 8pp. 

There the event is discussed in terms of the single production of a scalar top 
squark through a A' 131 coupling. The process e + p — > bxf — > bji+v^Xi with a sta- 
ble Xi would lead to a signature similar to that of our S6 channel. The process 
e + p — > bxt — » b^u^Xi ~~ > bfi + v^bdi> would lead to a signature similar to that of our 
S8 channel. 

[30] D.P Roy, Phys. Lett. B283 (1992) 270. 

[31] V. Barger, G. F. Giudice and T. Han, Phys. Rev. D40 (1989) 2987. 

[32] H. Dreiner and P. Morawitz, Nucl. Phys. B428 (1994) 31. 

[33] ALEPH Collaboration, D. Buskulic et al., Phys. Lett. B349 (1995) 238. 

[34] M. Drees and K. Grassie Z. Phys. C28 (1985) 451. 

[35] M. Hirsch, H.V. Klapdor-Kleingrothaus and S.G. Kovalenko, Phys. Rev. Lett. 75 
(1995) 17. 

[36] T. Kon, T. Kobayashi and S. Kitamura, Phys. Lett. B333 (1994) 263. 

[37] OPAL Collaboration, P. D. Acton et al, Phys. Lett. B337 (1994) 207. 



31 



